Buckling Instabilities of a Confined Colloid Crystal Layer 



0^ 



T. Chou and David R. Nelson 
Dept. of Physics, Harvard University, Cambridge, MA 02138 
(February 1, 2008) 

Abstract 

A model predicting the structure of repulsive, spherically symmetric, 
monodisperse particles confined between two walls is presented. When plate 
separations are small, only one layer of particles can be confined; however, 

> 

CN ' when the plate separation is increased, multiple layers will eventually form. 

o : 

. We study the buckling transition of a single flat layer as the double layer state 

' develops. Experimental realizations of this model are suspensions of stabilized 

cn : 

On ■ colloidal particles squeezed between glass plates. By expanding the thermody- 

' namic potential about a flat state of N confined colloidal particles, we derive 

a free energy as a functional of in-plane and out-of-plane displacements. As 
jH ■ the gap separation increases, certain out-of-plane modes soften. The wavevec- 

3 , tors of these first buckling instabilities correspond to three different ordered 



X 
J3 



structures. Landau theory predicts that the symmetry of these phases allows 
for second order phase transitions. This possibility exists even in the pres- 
ence of gravity or plate asymmetry. These transitions lead to critical behavior 
and phases with the symmetry of the three-state and four-state Potts models, 
the X-Y model with 6-fold anisotropy, and the Heisenberg model with cubic 
interactions. Experimental detection of these structures is discussed. 

PACS numbers: 82.70Dd, 46.30.Lx, 63.75.-hz 
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I. INTRODUCTION 



Colloid suspensions have long been used in practical applications such as paints, coatings 
and many manufacturing processes [0]. Recently, new uses of colloids such as in electro- 
optic devices and in semiconductors have been discovered 0. These new uses, as well as 
the general understanding of colloidal ordering applicable in other fields, require a thorough 
knowledge of the structures and phases present in these suspensions. Hence, colloid structure 
and rheology have been under intense study across many disciplines 0. 

Colloid particles, usually on the order of l/im in diameter, are a convenient size for 
laboratory study. Due to the complicated interactions these particles experience, colloid 
suspensions exhibit a rich cornucopia of phases in solution. These interactions can be de- 
scribed by electrostatic, van der Waals, hydrodynamic and steric effects. When repelling 
moieties such as polymers or charged molecules are attached to the surfaces, the colloids are 
stabilized against aggregation and extended ordered structures in solution can form. Debye- 
Huckel theory is typically used to model the screened electrostatic interactions in a charge 
stabilized suspension. Electrostatic effects, together with van der Waals and hard sphere 
interactions, are important ingredients in the Derjaguin-Landau-Verwey-Overbeek (DVLO) 
theories 

In bulk, monodisperse colloidal systems can appear in liquid, crystalline, amorphous, and 
inhomogeneous phases. Crystals with bcc, fee, and random stacked close packed structures 
have been observed. 

Theoretically, the ordering of hard or otherwise repulsive spheres have been approached 
analytically and via computer simulations. Much effort has been directed towards under- 
standing the nature and consequences of DVLO pair potentials. A Yukawa interaction was 
found to give rise to structural features consistent with observations . Pair potentials have 
also been used to described the dynamics of colloid crystals 0]. 

Experimentally, these systems have been studied with microscopy, video imaging, and 
light diffraction In bulk, monodisperse colloidal systems can appear as liquid, crys- 
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talline, amorphous, and inhomogeneous states. Crystals with bcc, fee, and random staeked 
elose paeked struetures have been observed [0,0- 

It was quickly realized that colloid suspensions are an ideal system to use for studying 
two dimensional phase transitions. The first observation of 2D ordering was the interfacial 
colloidal crystal of Pa. Pieranski |Q. In this system, colloid particles floating on an air- water 
interface interact mainly through dipole-dipole forces. Subsequent studies have involved 
immersed colloids with regions where two plates are brought close together. The plates are 
immersed in bulk solution allowing for free expansion of the confined structures. In addition, 
by using materials such as poly(styrene), which is almost density matched to the aqueous 
solvent, the effects of gravity are negligible. In these studies, the bulk solution provides a 
fixed chemical potential and thermodynamic particle resevoir for the confined sample. 

One class of experiments for which confined colloidal systems are particularly suited is the 
study of melting in two dimensions 0. Because of their macroscopic size, colloid spheres do 
not feel the atomic scale granularity of the confinement walls, (see Fig 1). Thus, this system 
is an ideal candidate for observing dislocation and disclination driven melting mechanisms 
in 2D. In fact, the search for defect mediated melting appears to be most successful in 
these confined colloidal systems. An isotropic fluid exists at narrow wall separations; as h is 
increased, the one layer areal density increases. In practice, h is difficult to accurately control 
and a wedge geometry with gradually increasing plate separation is used. The correlations 
of particle positions along this wedge are consistent with a two stage dislocation/disclination 
mediated melting process. 

If h is increased further, the continual invasion of particles into the gap is allowed only 
by formation of multiple layers. The observations of many experimenters can be sum- 
marized by, 

^ fluid ^ lA ^ 2n ^ 2A ^ ■ ■ -NU NA ^ {N + !)□ (1.1) 

where the gap h is increases from left to right and A and □ represent triangular and square 
lattices within each of the layers. In the gap, packing constraints dictate whether the 
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layers are square or triangularly ordered [|T^]. The □ and A layers are approximately the 
(001) and (111) planes of bcc and fee lattices respectively. As these transitions proceed, 
Murray et. al. have observed that large fluctuations of the colloids perpendicular to the 
plane occur; such fluctuations are especially prevalent in the lA ^ 2n transition. As 
increases, the structure of the confined layers approaches that of a bulk solution. In this 
study, we use static colloidal particle interactions to model the free energy of the assembly 
as a function of particle displacements. Stability conditions are imposed and minimization 
based on the Landau theory of structural phase transitions [|ll|701 is used to predict the 
possible buckled structures in the lA — > 2n region. The free energy depends upon both the 
amplitude and phase of these "buckling waves" . 

In the next section, we derive from microscopic interactions the free energy of the con- 
fined colloid system. Stability of this model and the critical wave vectors are examined in 
Section III. The phase transitions into each of the possible states is related to continuum 
statistical models in Section IV. Structure factors and other experimental consequences of 
these structures are then briefly discussed in Section V. Finally, a model for particles trapped 
in a tubular pore is presented in the Appendix. 



II. MODEL FREE ENERGY 

In this section, we derive a free energy functional from the thermodynamic potential of a 
confined layer of particles immersed in a bulk particle resevoir. The functional is written in 
terms of particle displacements and cast into the Landau-Ginzburg form. The parameters 
in this expansion are defined by the microscopic interactions in the problem. 

For a confined collection of A^ particles mutually interacting via a spherically symmetric 
potential U{r), the free energy can be written as a sum over particles i and j, 

^ = E U{\ri - rj\) + V{U h) - liN (2.1) 

where 
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V{f, h) = V,{^-f) + V2i^ + /) + VeM) (2.2) 

is the total particle-plate potential. Vi and V2 are the individual plate-wall potentials, 
which we usually assume equivalent when the confining plates are identical. The term 
Vextif) describes the effects of external fields such as gravity. In that case Vext oc f. The 
z component of Tj is denoted by fi, the height of particle i measured from the minimum of 
V^i + V2. 

For a single layer of spheres squeezed between the plates just prior to buckling, the 
mean particle positions form a flat triangular lattice with lattice constant a. Although in 
two dimensions there is only quasi- long range translational order |T^, we consider a flnite 
sample that is locally crystalline. 

Assuming our particle sample is within a well ordered domain, the sum over all nearest 
neighbor interactions becomes, 

Qo _ SNUjai) + 3NU{a2) + . . . - fiN + V{0, h)N 

Ao ~ N{a?^/2) ^ ■ ^ 

where Ao = Na?\/?)/2 is the area occupied by the undistorted A^ particle lattice. The last 
term varies with position in the wedge geometry of 0: this experimental arrangement is 
equivalent to slow spatial variation in the chemical potential. The second and further terms 
on the right result from summing interparticle energies from next and further neighbor 
particles respectively; hence, ai = a, 02 = oa/S, etc. For ( |2.3| ) to be useful, the further 
neighbor interactions must be truncated at a suitable distance much smaller than the A^ 
particle system size. Upon minimizing Vto/Ao with respect to a, we flnd the equilibrium 
lattice spacing and the chemical potential associated with it. 

The lattice spacing a is uniquely determined by the chemical potential provided U (r) is 
smooth and has at most one minimum. However, further neighbor interactions and interac- 
tions with multiple minima could lead to phase transitions among flat structures of varying 
lattice constants a. Multiple minima are often predicted by the various contributions to 
interparticle potentials. We will assume that the repulsive part of the potential is strong 
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enough to prevent aggregation and that there are no in-plane lattice distortions near the 
buckling transition when / becomes nonzero. The necessary requirements on U{r) will be 
determined below. 

Upon expanding the total free energy about the flat phase f — and about an equilibrium 
spacing a — |xj — Xj| between nearest neighbors with the definition of the in- plane 
displacement field u(xj+j) — u(xj) -|- a" = Xj — Xj, we have 

^ = E E E ^( V|u(x) - u(x + a-) - |2 + |/(x) - /(x + s.f)\^ ) + ^ V(/(x)) - /.(a) TV 

X n i=l X 

(2.4) 

Here, J2x is a discrete sum over the sites of a triangular lattice and n indexes different 
coordination shells. The {a"} are the three vectors {i} spanning the triangular lattice and 
inclined at 120° to each other, each with length a„, as depicted in Fig. 2. For nearest 
neighbors, 

—{^i } = {^i} = {y, - -y, -—X - -y}. (2.5) 

In order to write the free energy as a sum over a fixed number of particles, we consider a 
homogeneous system and a bulk chemical potential which, together with V, determines the 
lattice spacing and fixes the density of particles. 



A. Symmetric particle-wall potentials 

Upon expanding fl = fl/A about the height /=0 and the in-plane displacement u = 0, 
and using //(a) calculated to the appropriate order of further neighbor interactions, the 
nearest neighbor free energy density becomes, 

Q = const. + ^ EE ■ A.u(x)|^ + ^ E E (|A.u(x)|2 + |A,/(x)|^) - 

X i=l X ^=1 

^ EE(e. • A.u(x))|AJ(x)r + ^ EE |A./(x)r + 

X i=l X ^=1 

j:j:\Aj{^)\'+...+v{f,h) (2.6) 



X 1=1 



with 



V{f, h) = Y: (^f (x) + ufi^) + vf{^) + . . J (2.7) 



and 



(2.8) 



the projected area of colloid particles to second order in u. In these and subsequent 
equations, objects with direction indices {i} are cyclic: Oj = 0^+3. The notation a^^ = 
and Kl 23 = -^1.2,3 corresponding to a nearest neighbor approximation is used. In going from 
( p.4|) to ( p.6|) , terms linear in u have been cancelled by fi{a)/A. For further neighbors, we 
use the notation A"0 = 0(x) — 0(x + a"). To include the effects of the n^^ nearest neighbor 
spheres at in-plane distance a„, additional terms such as J^i.x I /(x)p must be added 
to Eq. ( p.6| ). The constant term will be omitted in future equations. 

The parameters K^, K2, and are completely defined by the microscopic particle- 
particle interactions, 

- -§M (2.9) 

an 

Kt^ ^ ^(a^) - Ki""^ (2.10) 



Kt^^^-f^^M-Kt^ (2-11) 



The expression for V is a Taylor expansion of the total particle-plate potential. In the 
absence of an external field and using identical confining plates, the function V{f) is sym- 
metric and its expansion contains only even powers of /. The terms in equation ( |2.6| ) contain 
the contributions from competing interactions. For repulsive particle-particle interactions 
Kl < 0; consequently the energy density maximizes |Aj/(x)p throughout sample. However, 
when repulsion from the walls, V{h, /), is strong enough, buckling is prevented. 



B. Asymmetric Potentials 



Thus far, we have considered only symmetric particle wall potentials in the absence of 
external perturbations. In general however, an external field may couple linearly to the 
height function /, 



(2.12) 



These perturbations may result from the effects of gravity or plate material asymmetry. An 
externally applied electric field would couple in this manner but would also induce dipoles 
resulting in an additional repulsive force between the particles. The total wall-particle 
potential becomes. 



V{f,h) = '-^f + uf + vf + ... + Xf 



(2.13) 



In mean field theory, this linear term shifts the minimum of each particle from / = to 
/ where / is found by minimizing ( |2.13| ) with respect to /. Keeping terms only to fourth 
order in /, we find. 



(3m 



,1/3 



9A2 + r^/3u - 3A 



2/3 



2(3u)2/3 yg\2 + rV3n - 3A 
Re-expansion of (p.2|) about / yields [|14 



— r 



1/3 



(2.14) 



V{f, h) ^ const. + I -Tji^f + A?7^(x) + Mr7^(x) + 

■V- \ ^ 



(2.15) 



Here, tj = f — f, and the new expansion parameters are related to the ones in symmetric 
potential case by. 



r = r + 12m/^ 



(2.16) 



A = Auf. 



(2.17) 

When asymmetry is taken into account, ( [^.6D is modified by replacing / with 77, r by f and 
adding the term Xi]^. This cubic term changes the nature of the transitions. 
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III. STABILITY ANALYSIS 



We now examine the region of validity of the free energy ( |2.6| ). In particular, we determine 
when and if phonons in the u or / fields become unstable or "soft". This analysis is most 
convenient in Fourier space. We define the Fourier modes, 

u{x) ='^u(k)exp{ikx.) (3.1) 

k 

and 

f{^) = J2f{k)exp{tk^) (3.2) 

k 

where /(k) = /*(— k) and u(k) = u*(— k) since the displacements are real. When periodic 
boundary conditions apply, the sums are over wavevectors spaced 27r/L apart where L is the 
linear dimension of of the particle system. 
To quadratic order, the free energy becomes 

/(k)/(G - k) + ^ |: u(k)- (y: M„(k)^ ■u(G - k) 

(3.3) 



OAT A ^ 



In (|3.3|) , we keep only the pieces quadratic in u which are not surface terms as they do not 
affect the eigenvalues cij^(k) of this dynamical matrix. The G are reciprocal lattice vectors 
with |G| = A7i/V3a. 

In this paper, we consider buckling instabilities in out-of-plane modes with the modes in 
u remaining stable. For the in-plane phonons, stability requires that the eigenvalues of the 
dynamical matrix, X]nM„(k), be real and positive. In the nearest neighbor approximation, 



Mi(k) 



Di(k) + §-^'2(2 — cos k-a2 — cos k-aa) ^_K'2(cos k-a2 — cos k-as) 

^ ^i^'2(cosk-a2 — cos k-aa) Di(k) + ^(6 — 4cos k-ai — cos k-a2 — cos k-aa) ^ 

(3.4) 



where 
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Di(k.) = 2Ki{3 — cos k-ai — cos k-a2 — cos k-as) (3.5) 
We find that the all eigenvalues of Mi (k) are always real and are positive provided, 

K, + ^>0 (3.6) 

and 

1 1 \^ 1 ^ 

^ 2KI) ^ Df{k) ^ (^^^^ ~ ^^'^ k-a,+i) (3.7) 

with SLi = aj+3. The criterion ( |3.6D is satisfied with repulsive potentials in general. In the 
first Brillouin zone, the right hand side of ( ^.71 ) varies between and Ki^/16. Therefore in 
the nearest neighbor approximation, in-plane phonons are stable provided, 

K, + ^)(k, + '-^)>0 (3.8) 



4 / V 4 

This is just a more stringent version of ( |3.6| ) and satisfied by potentials which are repulsive 
and short ranged consistent with the nearest neighbor approximation. 

Now, we find instabilities in the magnitude of /(k) and assume that the phase of the 
displacement functions have been fixed to minimize the free energy. The dependence of the 
corrugations on the phase of / will be discussed later. Provided (|3.71 ) holds, we can consider 
only the instabilities in the transverse modes, /. Here, we require 



UJ 



■(k) = r(/i) + 5]Z}„(k)<0 (3.9) 



for instability of the k*'' mode of /(k). As the spacing between the plates increases, the 
compressional forces balancing the inward forces from the colloid bath and r{h) decrease. 
The points at which the criterion ( p.9|) is first satisfied selects the first modes to soften. An 
estimate of the critical spacing he can be made by solving r{hc) — — Z]n-D„(kc). 

In the nearest neighbor truncation, r{h) = — Di(k); the values of k^ that satisfy this 
relation for the largest value of r{h) are at the corners of a hexagonal Brillouin zone shown 
in Fig. 3(a). Thus, as r{h) is decreased the first out-of-plane corrugations of the colloidal 
array will occur at these wavevectors. When further neighbors are considered, kc can shift 
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to the edges of the Brillouin zone, as in Figs. 3(b). The location of these first instabihties 
depends on the relative sizes of k[^\ In the next nearest neighbor truncation, the first 
unstable modes occur at the zone edges when 

\K?^\>l\K,\ (3.10) 

o 

as shown in Fig. 3(c). With repulsive and reasonably short ranged potentials (such that 
^ 0), the corner and edge instabilities are always the ones to occur first. We will 
restrict ourselves to these two types of instabilities. The real space structures of these 
corrugations are depicted in Fig. 4. 

Since three pairs of corner instability points are connected via reciprocal lattice vectors, 
the height function / is the real part of a one component complex function. The edge modes 
can be represented by a three component function. The sums over critical wavevectors kc 
are thus restricted to the stars of ±ko and q^, {i = 1,2,3) belonging to the symmetries of 
the corner and edge instabilities respectively, 

An 

ko = ±—y (3.11) 
3a 





27r ^ 




qi = 


\/3a 










q.2 = 




-y 

a 








qs = 


—TT ^ ^ 


-y 

a 



(3.12) 

These two sets of unstable modes transform as particular irreducible representations of 
the space group of the lattice. Since the height function can be written as linear combinations 
of the basis functions of one of these two irreducible representations (excluding the unit 
representation), a continuous transition into these ordered structures is possible if no cubic 



invariant in the free energy exists |TT| , p!5| . This Landau rule is not exact in two dimensions 
because fluctuations can lead to continuous transitions even when the cubic term is present 
I^,|2D[|. Landau theory nevertheless plays an essential role in determining the universality 



class of the transition. The functions /(±ko) and /(q^) are identified with d component 
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order parameters where d = 2,3 are the dimensionahties of the irreducible representations 



15,16 



When unstable modes with wavevectors at zone boundaries occur, the energy will depend 
on the phase of /. For example, with k and — k pointing to the midpoint of an edge of the 
Brillouin zone boundary, careful enumeration of the sums over wavevectors of the reduced 
group, or star, leads to a quadratic term oc \{fkfk + c.c). 

The free energy without phase terms is correct for wavevectors such that k + k' 7^ G 
and is sometimes termed the "incommensurate" free energy. However, when the relevant 
wavevector is is a rational fraction of the reciprocal lattice vector, the phase variable enters 
and a "commensurate" free energy obtains |jl2|]. Since the functions I]n-^n(k) for n < 3 
appear to have stable extrema only at corners or midpoints of edges, the commensurate 
wavevectors (p.ll|) and ( |3.12| ) are the most important. Thus we expect the system to order 



directly into a corrugation with commensurate wavevector. 

We now derive these commensurate energies. Since the energy was a sum over discrete 
lattice sites, the expression (^.31) contains sums over G. Thus, at the boundaries of the 



Brillouin zone, certain wavevectors can participate in the "Umklapp processes" shown in 
Fig. 5. It is these terms that generate dependent terms in the energy that lock the 
phase. The phase dependence can be incorporated by generalizing the order parameter, and 
allowing for a slow spatial undulations in both the phase and amplitude, 

/(k) = Mx)e''f'-^^\ (3.13) 

where f\i{x) and 0k (a;) are real positive functions. The procedure we follow is first fixing 
the amplitude and phase of the order parameters, expressing the energy in terms of these 
variables, then allowing slow spatial modulations within the sample. 

At kc = ±ko (the corner instabilities shown in Fig. 3(a)), the nearest neighbor, u(k) 
independent part of the commensurate free energy becomes, 

{r{h) + 9K,) \f{K)\' + \(Qu+ ^) \f{K)\' + +j (2OV + |/(k,)r 
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+1 (2^ - ^) co.(60)l/(ko)r + • • • (3.14) 

Since 2v — Slf^a/Sa^ > 0, the six phases minimizing (|3ll) for the 73 X V3 state are 
{0m} = {(2m + l)vr/6, m = 0, ...,5}. They represent three equivalent states generated by 
translations in three directions. This degeneracy in each of these states is doubled by the 
inversion /(k^) -/(k^). 

For edge wavevectors the quadratic terms fix the phase of the displacement fields at 
= mvr, provided higher order terms are sufficient to stabilize the energy. For these critical 
modes, kc = (i = 1, 2, 3), we have 



j + ^ + ^ ) E [3i/.n/mpco.2(0, + 0,+o + mn\f^-l\' + i/mnco.20,] + . . . (3 



where we have set fi = /(qi) and 0j = 0(qi). In this case the quadratic term selects the 
phase. For AKi + AK[^^ + r/2 < 0, = mvr and the quadratic coefficient in ( p.l5| ) at this 
phase is 8Ki + 8K[^^ +r{h). Comparing ( p.l4|) with (|3.15|) leads to the condition ( |3.12|) for 
edge instabilities. 

When edge instabilities and phases lock as described above, the order parameter still 
contains three independent real components, {|/i|} corresponding to three equivalent direc- 
tions for the corrugation. In principle, the quadratic couplings in Eq. ( |3.15| ) could select 
a "triple-q" structure (see Fig. 4(c)), in which all three components order simultaneously. 
However, in the symmetric wall potential case, the ordered state has a configuration with all 
the weight in one of these components which is spontaneously chosen by the system. This 
"single-q" state minimizes the free energy for symmetric wall potentials and is depicted in 
Fig. 4(b). This resembles an observed maze-like pattern |p|,p!7| and can be an intermediate 
structure in the transition to the 2n phase. However, we do not expect antiferromagnetic 
Ising behavior (see below) as suggested by Ogawa [|17|. For the free energy (p.l5|) , a 3-q 
state, where all the \fi\ are equal, is favored only in the presence of asymmetry in V{f) 
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IV. CORRESPONDENCE OF FREE ENERGY TO CONTINUUM STATISTICAL 

MODELS 



We now explicitly allow the order parameter to acquire a slowly spatially varying modula- 
tion, as defined in ( |3.13| ). This is equivalent to expanding k-dependent quadratic coefficients, 



Di(k) about k ~ ko. With the notation fo{x) = fkoi^) and (po{x) = 4>]f_^{x), the nearest 
neighbor free energy in position space becomes, 

3a2|i^,| , 1 , ^,,,,,2,., , 1 , 81 



^ = ^J^l V/o(x)|^ + - {9K, + r{h)) f^{x) + j[Qu + —K,j f^{x) + 

J (2OV + ^i^a) foix) + ^^3) cos 60,(x) + ■■■ (4.1) 

At a fixed amplitude which minimizes fo, the expression describing the energetics of the 
phase variable becomes 

- = ^^|V0,(x)|Vo' + i - foCosGMx) + ... (4.2) 

This expression has the form of the 2D X-Y model with six-fold anisotropy as discussed 
by Jose and others p3|JT8[| . This model is expected to order in two steps with an intermediate 
region of continuously varying exponents as the parameters are varied. 

When asymmetry exists, 

^ = ^4|^|V0o(x)|\2 + jXvlcosSMx) + ... (4.3) 

Here, the largest phase locking term is proportional to rj^cos 30o To minimize the energy, 
the phases prefer the values 0o = {if?^} when A > and {0,±^} when A < 0. This 
hamiltonian is identical to that of the continuum 3-state Potts model [|1^]. The behavior of 



this hamiltonian has been studied extensively p!9|-p2|; a continuous transition is possible, 
and the critical exponents are known . 



When the edge modes buckle, there exist two cases. If a 1-q state is favored, any one of 
three directions is selected, each with two preferred phases 0. These conditions are described 
by a Heisenberg model with quartic couplings proportional to {j2i=i l/iP) and J2i=i \ fi\^- 
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The exact critical behavior of this model in two dimensions is, to the best of our knowledge, 
as yet unknown. 

In fact, for the particular hamiltonian at hand, the 3-q state is energetically preferred 
over the 1-q only when appreciable symmetry breaking interactions (gravity or asymmetric 
plates) are present. The term, 

3 

A J2 Vimm cos(0i + - 0i+2) (4.4) 

i=l 

favors the 3-q or 2 x 2 state shown in Fig. 4(c); energy configurations with any rji = 
will be unaffected. In the presence of such a term, the degeneracy in phases of the 

3- q structure is halved. Now the sets {0j} = {(vr, vr, vr), (0, 0, vr), (0, vr, 0), (vr, 0, 0)} and 
{(0, 0, 0), (tt, TT, 0), (tt, 0, tt), (0, TT, tt)} are selected when A > and A < respectively. A 

4- state Potts model (which also has a continuous transition in d = 2) is recovered in this 
case [^,^,^. Of course, as r{h) is decreased further, discontinuous transitions from the 
3-q to 1-q corrugations are possible. 

The minimum energy configurations in order parameter space are shown in Fig. 6. In 
general, a state described by an equally weighted superposition of directions, the 3-q state, 
has eight possible combinations of the three phases associated with each order parameter 
component that minimize the energy, as shown in Fig. 6(c). The -\/3x ^/3 corner instability is 
represented by Fig. 6(a) and Fig. 7(a) for the symmetric and asymmetric cases respectively. 
For the free energies derived from our microscopic model, the situation depicted in Fig. 6(c) 
does not occur since in the absence of a symmetry breaking term, the 1-q state (Fig. 6(b)) 
is always preferred. 

V. EXPERIMENTAL CONSEQUENCES 

In this section, we briefly examine the experimental signatures expected for each of the 
buckled structures above. To begin, we must have a rough description of the interaction 
parameters. The simplest DVLO type potentials for the sphere-sphere and sphere-wall 
interactions in the charge stabilized limit are the Yukawa and exponential forms respectively, 
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U(r) = —e-^' (5.1) 
r 



In the experiments of Murray et.al. 0, a layer of O.Sfim poly(styrene) sulfonate spheres 
were confined in a gap h ~ 3^m. To stabilize the structures, low ionic strength solvents 
with K,a ~ 0.05, where n is the Debye screening length, were used. With these estimates we 
find \Ki\ ~ 1 X 10~^erg/cm'^, r{h) ~ OAerg/cm^, and u{h) ~ lerg/crn^. Interactions of the 
form ( p.2| ) make r and u very sensitive to the gap spacing h and very small wedge angles are 
probably needed to study these transitions in detail. Although our analysis has neglected 
the equilibrium thermal fluctuations of the spheres in solution, it will correctly identify the 
symmetry or "universality class" of the various transitions. Thermal fluctuations will change 
the detailed predictions of mean field theory whenever the mean square fluctuations of the 
order parameter amplitude exceeds its equilibrium value. 

Asymmetry effects can be described by the additional term ( |2.12|) . If the external field 



of concern is gravity A = fo(pc — Po)g, where Vo is the volume of each particle and pc and po 
are the densities of the colloid and suspending fluid respectively. For poly(styrene) spheres, 
A ~ 7 X lO^^'^ erg /cm. 

An estimate of the relative displacement of the entire layer from the center of the gap 
using the values in (|2.14D gives \ f /h\ ~ 1 x 10~^ which is truly negligible. Another measure 



of the degree of asymmetry is the change in the critical value r*{h) where the magnitude /k 
first becomes non-zero. At a corner instability (-\/3 x -\/3)j for A << Kf/u, 

Ao?u 



4A2 

r* — Sr* ~ 



24ua^ + 8IK2 



(5.3) 



Since K2 > and / depends on r{h), the critical value r*{h) is lowered by the presence of 
A, i.e., a larger h* is required to buckle the lattice. 

Given that the conditions for equilibrium ordering are met, there are essentially two 
methods currently used to probe these structures, direct video imaging and diffraction. 
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When thermal fluctuations are strong, laser light diffraction patterns can be analyzed to 
elucidate average structural features. 

For monochromatic light incident on the layers (the plates are made transparent), diffrac- 
tion peaks are expected at wavevector transfers equal to the critical wavevectors ( p.ll| , prT2D . 



The 73x^3 structure is defined by a primitive cell with three times the area of the original, 
undistorted lattice. A three point basis is also superimposed. The static, zero temperature 
structure factor for point, single scatterers in this case is, 

5(kx, G)5(-kx, G') [2 + 2 cos{k ■ ei + hzo) + 2 cos(k ■ (ei + 63) + hzo)] (5.4) 

where the 2D reciprocal lattice corresponding to the larger hexagonal lattice are, 

r^T f nAir ^ mAn.y/S^ 1 , ^. 

{G} = ^y, ^h^x + -y 5.5 

and Zo is the spacing between the two flat lattices, Zq = |/kc- In a x \/3 structure, 
in addition to the six spots expected from the original triangular lattice, we would see six 
additional first order peaks at smaller wavevector transfers and rotated by 30" provided Zq is 
large enough. This pattern is in marked contrast to the asymmetric, three spot pattern seen 
in diffraction from a 2 A superlattice. We note that all diffraction spots in two dimensions 
are algebraic singularities |jl3|, rather than true (5-function Bragg peaks. 

The 2x2 structure has even a larger primitive cell and a four point basis. The six 
diffraction pattern peaks would be even closer to the origin that in the x ^/3 structure 
and would not be rotated with respect to the six original peaks. A pattern consisting of 
four spots rectangularly positioned would be an indication of a 2 x 2 superlattice. A similar 
square pattern develops when the ordering is 2n. It would be interesting see how these 
patterns evolve as r{h) is decreased. 

In addition to diffraction patterns, video imaging of the colloidal crystal may reveal 
the out-of-plane buckling and the concomitant in-plane lattice distortions. The in-plane 
contractions are due to the ei ■ Aju| Aj/p terms in which we have neglected. To predict 
any effects in the lattice spacing as buckling occurs, we minimize the energy density, Q with 
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respect to the AjU associated with the specific corrugations in Fig. 4. The AjU used are 
the k = modes expected on symmetry grounds for each of the structures. A x -\/3 
buckhng is expected to induce a symmetric contraction of magnitude, 

The contractions depend on (/^) which, in turn behaves according to the relevant model 
near the transition with a |r — r*|^~" energy-like singularity. Similar expressions hold for the 
other buckled structures. Direct observation of lattice contractions may reveal the nature of 
the buckling transition even if the predicted extra spots are blurred by thermal fiuctuations. 
It is not clear how strongly the in-plane displacements are coupled to the buckling at larger 
amplitudes. For instance, individual fiuctuations in /(x) of the charge stabilized particles 
can induce a fiuctuating dipole in the z direction which would contribute a net repulsion 
between the particles tending to expand the lattice. 

VI. SUMMARY 

Stability analysis has shown that the buckling of a layer of confined repelling particles 
results from the competition between a bulk compressional force derived from a chemical 
potential and a confining wall force derived from a repulsive particle-wall potential. We 
have found that the unstable modes have wavevectors kc at the midpoint of edges and at 
the corners of the Brillouin zone boundaries. These instabilities correspond to the 2 x 1 or 
2x2 and the x structures respectively. In this analysis, the edge instabilities arose 
only when further neighbor interactions were considered. The additional restriction ( p.9|) 
was imposed for the in-plane phonons to remain stable. 

The symmetry allowed buckled states have basis functions belonging to irreducible repre- 
sentations of the point group. Furthermore, an expansion of the order parameter coefficients 
in the wavevector index k yields minima at kc. Despite the Landau rule prohibiting a con- 
tinuous transition in the presence of gravity or plate asymmetry [|r^,|r^,|T^] , we find that in 



this two dimensional problem, second order transitions may occur [p!q , P^ . 



Provided we allow the order parameter phase to be slowly space dependent, the free 
energy can be written in a continuum form. A striking feature is that these hamiltonians 
allow for continuous transitions into the mentioned buckled structures, and are predicted 
to have critical behavior corresponding to a X-Y model with six-fold anisotropy for corner 
instabilities, and a Heisenberg model for an edge instability. Colloids between symmetrical 
plates may be the first experimentally realized nonmagnetic system which can be described 
by an X-Y model with six-fold anisotropy. In general, the edge instability has critical 
behavior resembling that of a Heisenberg model with cubic interactions, although in this 
particular problem transitions to 3-q states are expected to occur only in the presence of 
asymmetry. 

Asymmetry from the effects of gravity were found to be negligible. When asymmetry 
is expected to be important continuous transitions are not precluded by cubic invariants 
in the free energy. In these cases the hamiltonians are in the 3-state and 4-state Potts 
universality classes for the corner and edge instabilities respectively. It would be interesting 
to use differently treated upper and lower plates to enhance this asymmetry. 

The structures enumerated above should be easily distinguished using light diffraction. 
We have estimated the structure factors and predicted the diffraction patterns that would 
result from these superlattices. It is also possible to obtain an estimate of the contrac- 
tion in lattice spacing due to the buckling. This effect is described by the coupling term 
(e,A,u)|A,/|2, inEq. 

The main prediction of our theory lies in the analysis of the types of instabilities that 
may occur. We emphasize in conclusion that the precise mode of buckling and hence which 
type of phase transition occurs depend sensitively on the microscopic parameters in the 
theory. 
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APPENDIX A: ONE-DIMENSIONAL SYSTEM 



In this section, we consider a related problem: a one dimensional string of confined 
particles as shown in Fig. 8 For a line of regularly spaced nearest neighbor repelling particles 
confined to a cylindrical tube, the free energy expression equivalent to (|2.6| ) is. 



n 



[Ki + K2) 
2L 



N 



K ^ K 



N 



1 J^/r, 



(Al) 



where the length of particles is exactly. 



Na 



N 



1 + E(^i+1 ~ 



(A2) 



Equation ( [Al| ) can be supplemented with similar terms to include further neighbor inter- 
actions. The parameters r{d) and u{d) may be tuned by changing the pore diameter, or 
alternatively, by using a cone geometry. In contrast to the two dimensional sheet, the in-line 
displacements Ui are scalars and the out-of-line displacements fi are two component vectors. 
The harmonic free energy expression analogous to ( p.3|) is. 



Mn{k) = + irr)(2 - 2cosnka) 



(A3) 



and, 



{d) + Y,K^r\2-2cosnkc 



where i,j label two orthogonal directions perpendicular to the tube axis. 



(A4) 
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The flat state in the nearest neighbor case ( |A1|) is unstable at wavevectors kc = ivr/a. 
The free energy of a structure with an instabhlity at /c^ = vr/a becomes, 

" = + + - ^'^*'') ^" ^ ■ ■ ■ 

In one dimension, however, the inclusion of the next nearest neighbor interactions pushes 
the first unstable modes inside the Brillouin zone, to smaller k^. The free energy of these 
longer wavelength modes will have higher order phase locking terms only if kca/n is rational. 
Otherwise, the free energy resembles that of an isotropic X-Y model, 

^ = ^ E Dn{kMkc)\' + ^ E E ^^«V^(M^^|^|V0.(x)r + . . . (A6) 

kc,n kc,ni=x,y 

Furthermore, each component of fk contains a phase variable. In addition to the absolute 
phase, the phase difference, \(f)x — 4>y\ also determines the structure. At zero temperature. 



the phases minimizing the nearest neighbor approximation ( Al ) are 0j = mvr and the phase 
differences are or tt corresponding to a linearly polarized wave as shown in Fig. 8a. 

If, from the inclusion of further neighbors, kc is incommensurate, or if the first phase 
locking term is very high order, the prefered phase differences \(j)x — vary continuously. 
For example, when 4:\Ki \ > \Ki\ the next nearest neighbor truncation, (n = 1,2), yields, 

cos krtt = 7^ (A7) 

and, 

(r{d) + 2K, + AKf^ + + E «Vf (A:.) (-^y- - ^) |V0.(x)p + 

^ \ 4:Ki V i=x,y \2Ki 'L ^ J 

(A8) 

The mean field structures available range from linearly polarized undulations {{(^^ — (pyl = 
mvr) to circularly polarized helices {\(f)x — 4>y \ = (2m + l)/27r). At low enough temperatures, 
these spontaneously chosen motifs may be observable. 
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FIGURES 



FIG. 1. A confined layer of repelling colloid particles. 

FIG. 2. Vectors indexing coordination shells for a triangular lattice, la^-*^)] = a; |a(^)| = ^/Sa, 
|a(2)| = V3a,... 

FIG. 3. The function w^(k) = r{h) + Di(k.) + D2{'k) slightly above and exactly at the transition. 
The first Brillouin zone is shown in the inset. Zeros first occur at a corner instability (a) when 
\k[^^\ > ^\Ki\ or (b) edge instability otherwise. 

FIG. 4. Minimum energy structures with out-of-plane displacements indicated, (a) -\/3 x -\/3 
corrugation, (b) 2 x 1 structure, (c) 2 x 2 structure. Note that not all particles have equal 
displacement magnitudes; however, /(x) = 0. 

FIG. 5. "Umklapp" processes for corner and edge instabilities. In the absence of asymmetry, 
G = 2qj for the edge modes: the phases are fixed by the first quadratic term. For the corner 
instabilities, G' = 6ko is the first Umklapp term. When a cubic term exists, a lower order Umklapp, 
3ko = G" appears. 

FIG. 6. Order parameter space representation of model Hamiltonians. Energy minima are 
shown by filled circles, (a) The X-Y model with six-fold anisotropy when corner instabilities 
arise, (b) Heisenberg model when a 1-q state is chosen within the edge instability subspace. (c) 
Heisenberg model when "cubic" interactions dominate, (3-q). 

FIG. 7. Order parameter representation with asymmetry, (a) 3-state Potts, (b) Heisenberg 
model with cubic interactions. The light and dark points refer to opposite signs of the asymmetry 
parameter A. 

FIG. 8. A confining geometry of particles in a conical pore 
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